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ABSTRACT 

A review of some of the current literature pertaining 
to this thesis is conducted. The Created Response Surface 
Technique of solving a class of nonlinear mathematical 
programs is presented. The theoretical interpretations 
of the primal and dual formulations of the technique in 
an activity analysis context are developed. The applica- 
bility of these interpretations to the neoclassical theory 
of the firm and the contemporary organization theory is 
indicated. Computational experience in solving well de- 
fined numerical problems is also indicated. Several lin- 
ear models of organizational decentralized planning are 
vi # a mo d.0 i cf w0C0D."tT , 2.1.d.z0ci pic. 

procedures is proposed and solved using the Created Re- 
sponse Surface Technique. 
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I. INTRODUCTION 



, A. REVIEW OF NEOCLASSICAL THEORY OF THE FIRM 

In this thesis the neoclassical theory of the firm and 
contemporary organization theory are extended by the appli- 
cation of a nonlinear programming technique known as the 
Created Response Surface Technique (CRST) . For ease of 
understanding, the reader is expected to have a basic 
grasp of the neoclassical theory of the firm, the static 
economic theory of market equilibrium, and from the mathe- 
matical disciplines the classical theory of optimization 
including Lagrange multipliers and mathematical program- 
ming. This introduction contains a brief review of the 
neoclassical theory of the firm and a description of the 
background motivating this thesis. 

In neoclassical economic theory a firm, or producer, 
is an economic agent who consumes some commodities and pro- 
duces other commodities. It is usual to embed the pro- 
ducer in an economy characterized by scarce resources. By 
the interaction of consumers and producers in markets, 
commodities which are goods or services are valued as 
scarce, free, or noxious. More specifically, the firm 
technologically transforms inputs, factors of production 
which are purchased in markets, into outputs, products 
which are sold in markets. In making choices the firm con- 
siders consumers' and other producers' demand for its out- 
puts, the state of production technology and the behavior 
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of the market where it purchases its inputs and sells its 
outputs. The firm is usually assumed to behave as if it 
maximized its profit. In this case profit is the differ- 
ence between the total revenue received from the sale of 
its outputs and the total costs incurred in the production 
of those outputs. 

In a perfectly competitive economy a market mechanism 
establishes prices for all commodities including the firms 
inputs and outputs. The firm is assumed to accept these 
prices as given datum determined in the market (atomistic 
competition) and to pursue some form of optimizing behav- 
ior. This behavior can be described by three different 
behavioral models. First, there is the maximization of 
physical output subject to a cost constraint. Second, 
there is the economic dual formulation of that model, the 
minimization of total cost subject to an output constraint. 
And finally, there is the maximization of profit, which 
under suitable conditions embeds the previous two formula- 
tions as subproblems. 1 

If the production, cost and profit functions are suf- 
ficiently well behaved, 2 mathematical programming techniques 
can be applied to these models. The application of mathe- 
matical programming techniques to these problems usually 



1 

Henderson, J. M. , and Quandt, R. P., Microeconomic 
Theory , p. 53, McGraw-Hill, 1958. 

2 

Kuhn-Tucker Constraint Qualification, Reference, 
Appendix A, paragraph 9. 
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results in decision rules which characterize an optimum. 

If the functions or their approximations are sufficiently 
mathematically tractable, numerical solutions can be 
achieved. 

B. DECENTRALIZATION PLANNING PROCEDURES 
1 . General Model 

In neoclassical theory the firm is modeled as if 
it were an homogeneous entity. Mathematical functions 
represent the firm's aggregate behavior ex cathedra in the 
economy. The decision rules derived from these models in- 
dicate the necessary and sufficient conditions to achieve 
optimal behavior. 

Several models have been proposed which describe 
how the firm pursues the optimum behavior once the neces- 
sary and sufficient conditions for that optimum behavior 
have been established. One model proposed by E. Malinvaud 
[1]* considers the firm's pursuit of optimum behavior to 
be a sequential decentralized planning procedure. This 
model is representative of a general class of decentral- 
ization models and will be presented to introduce the 
basic concepts and the pertinent assumptions underlying 
such decentralized planning models. 

This model characterizes a firm which is suffi- 
ciently large to have such organizational subunits as a 
central agency and at least one production agency, denoted 



* References in brackets may be found in the List of 
References . 
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It is supposed that the central 



by j where j = 
agency is responsible for formulating a plan for opti- 
mizing the behavior of the entire firm. Assume that the 
central agency can forecast the aggregate demand for the 
firm's products over a time period in which technology 
may be considered constant. Further suppose that the cen- 
tral agency knows the quantities of available inputs and 
the entrepreneur's preferences, which may be represented 
by a utility function, but it does not know the specific 
technical information which describes the actual produc- 
tion process. This specific information is known by pro- 
duction agencies, but each specific production agency does 
not know the entrepreneur's preferences, the availability 
of inputs , or tlio pi*o duct ion isvols of otlion production 
agencies, if they exist in the firm. Clearly, some method 
should be used to allow the production agency to provide 
its specific technological information in the central 
agency's formulation of the plan. Although the formula- 
tion of the plan is in reality a continuous process, it 
is assumed that this plan will be created over a time 
period in which these exchanges of information between 
the central agency and the production agencies will form 
an iterative procedure . 3 

This iterative procedure is assumed to be accom- 
Pl ished in discrete stages. Each stage is denoted by s, 

3 

This conclusion presupposes some form of administra 
tive regulation governing how information is exchanged and 
some requirement for minimum content of the information ex 
changed . 
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s = 1, . . . ,S. The central agency provides "prospective in- 
dices", I s , a vector consisting of J components at stage s 
Each component of this vector is the "prospective index" 
to each of the J production agencies at the s^* 1 stage. A 
prospective index is essentially a request for technolog- 
ical information. It's form may be a production quota, 
resource availability, or a budget constraint. For ex- 
ample, the central agency may be asking, "What can you, 
the j th production agency, produce if you have "X" dollars 
of budget available?" Each prospective index requires 
some response from each of the production agencies. In 
return, the production agencies provide a "proposal", 

P s , a vector consisting of J components at stage s. Each 



component of this vector i 
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J production agencies. A proposal implicitly provides the 
specific technical information embedded in the produc- 

tion agency. For example, the j*-* 1 production agency may 
be responding, "Given "X" dollars of budget, I can produce 
"Y" units of my specific output." Each mutual exchange of 
information is a stage in the formulation of the "plan". 

The vector, p s , will denote the "plan" after the s" 1 ^ stage, 
This procedure is represented as follows: 



Central 
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where , 








2 . Characteristics 

Analogous to the characteristics of the simplex 
algorithm of linear programming, any iterative procedure 
can be studied to understand which characteristics it ex- 
hibits. The folloitfing is a partial list of possible 



characteristics 
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cedure could be well defined. That is, if at each stage 
of the process solutions exist to the "prospective indices", 
"proposals" and the "plan", the procedure is well defined. 
Secondly, the procedure could be monotonic. If the utility 
of the s^ iteration is at least as great as the utility of 
the s-ist iteration, the process will be monotonic. Next, 
the procedure could be convergent. If, as S increases 
beyond bound, the utility of the s^ iteration approaches 
some least upper bound of the set of utilities defined over 
all feasible plans, the procedure will be convergent. Fin- 
ally, the procedure could be finite. If the "plan" after 
the s th iteration is the optimal plan for some finite num- 
ber S, the procedure will be finite. If the procedure is 
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not finite, it is assumed that the central agency will 
establish an upper bound on the number of iterations of 
the process. In reality, the time and cost of each ex- 
change of information will tend to constrain the process. 
Some optimization technique is indicated and would be 
appropriate. This optimization technique will not be 
discussed, but it may provide an area for future interest. 

3 . Walrasian Tatonnement Process 

As an example of iterative planning procedure, 
consider the neoclassical Tatonnement (recontracting) 
process of static market equilibrium. A perfectly com- 
petitive market might consist of an auctioneer, customers, 
and producers (production agencies) who are assumed to 



produce only one uniq 
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the consumers and producers yields market equilibrium (the 
plan) by the Walrasian price adjustment mechanism. The 
price adjustment mechanism acts as * a stopping rule signi- 
fying the termination of the iterative procedure and 
achievement of the necessary and sufficient conditions for 
market equilibrium. The procedure starts when the auc- 
tioneer, real or pedagogical, provides a price vector 
(prospective index) of all the commodities in the market 
to consumers and producers. This initial vector may be 
based on past performance of the market or "judgment". 

The consumers seek to maximize their individual 
utility subject *to a budget constraint generated by the 
prices of commodities in their prospective commodity 
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bundles and their initial endowment. The quantities of 
goods and services the consumers are willing to buy at 
the existing prices form the aggregate demand at this 
iteration . 

Concurrently, producers seek to maximize their 
individual profits. Since profit is a function of inputs 
and outputs, the existing prices allow firms to determine 
the quantities of their respective outputs they are willing 
to produce. The vector of these quantities (proposal) is 
the aggregate supply at this iteration. 

If consumer demand for the j t * 1 product at the s^ 

iteration, denoted by D . , exceeds production supply of 

the jth product of the s^h iteration, denoted by S ., or 

^ 3 

vice vers 3 . foT 3.ny of "tho commocli "txos , th.6 ductionooT* v/i.11 
adjust the appropriate prices to account for the excess 
demand or supply. The iterative procedure continues until 
the market is cleared, D s - = (T, at an equilibrium 

price, indicating that the plan is optimal. 
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II. REVIEW OF LITERATURE 



Activity analysis is the synthesis of mathematical 
programming and the economic theory of the firm. When 
mathematical programming techniques are applied to well 
defined behavioral models of the firm, the optimum solu- 
tion to the programming formulation of that model yields 
appropriate managerial decision rules which characterize 
that optimum solution. Likewise, if numerical data is 
used in a practical model of the firm, the optimum solu- 
tion to the programming formulation of that model will 
indicate how that optimum solution may be achieved with 
that data. 

Historically the theoretical and practical applica- 
tions of activity analysis to well defined resource al- 
location problems appear to have developed as a hybrid of 
research in both disciplines. The introduction of a 
mathematical programming technique is generally followed 
by the application and interpretation of that technique 
in an economic context. 

Whether or not a cause-effect relationship actually 
exists is not conjectured. Rather, this apparent sequence 
will be used to provide topic continuity in reviewing some 
of the principal research efforts in areas which directly 
relate to this thesis. 

The first known model of an economy which introduced 
technological production in characterizing a general economic 
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equilibrium was developed by J. von Neumann [2] in 1937. 
The model treated technological coefficients, outputs, and 
production processes as known entities. Activity levels, 
prices, the coefficient of expansion of the entire economy 
(scale of the economy) and the interest factor of the 
economy were considered unknown entities. The model could 
accomodate the concepts of joint production, capital goods 
depreciation, and scarce or free goods. Employing the 
Minimax criterion of Games Theory, von Neumann proved the 
existence of at least one solution to the coefficient of 
expansion and the interest factor as functions of prices 
and activities. At optimality, the solution to the model 
uniquely equates the coefficient of expansion to the in- 
terest factor. The equality of these implicit functions 
of prices and activity levels characterize the equilibrium 
conditions of the economy in perfect competition. 

During World War II the allocation of resources in 
multi -theatre military operations and in the domestic 
economy generated interest in mathematical models of de- 
finitive planning procedures. A landmark in the research 
effort in linear function analysis was G. B. Dantzig's [3] 
introduction of the simplex algorithm in 1947. 

T. C. Koopman's [4] immediate appreciation for the 
variety of allocation problems which could be modeled by 
procedures similar to the simplex algorithm resulted in 
the emphasis on activity analysis at the Cowles Commission 
for Research in Economics Conference of 1951. The primary 
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purpose of this conference was to consolidate the results 
of various independent research efforts which had been 
addressing similar problems in an activity analysis con- 
text . 

In a recent book, D. C. Vandermeulen [5] describes 
the foundations and origin of the simplex algorithm and 
its application to the neoclassical theory of the firm. 

The iterative procedure of the algorithm is interpreted 
as a model of the rational economic behavior of an hypo- 
thetical entrepreneur. The decision rules governing 
optimization of a defined production objective are dem- 
onstrated. The economic analysis of complementary 
slackness conditions and the applicability of the dual 

■Tn vm ill n -f* n a -p a m a a 1 a a 

xv jliuUx (i^xwu wo. wiiu iiiwuu j. a j. ^ ^ico^illuUt 

In 1950 H. Kuhn and A. Tucker [6] formulated the 
necessary and sufficient conditions which characterize 
the saddle point of a two person zero sum game and applied 
those conditions to a class of nonlinear inequality con- 
strained programming problems. Under appropriate dif- 
ferentiability and nonnegativity assumptions, it was shown 
that the maximum of a concave objective function over a 
convex set defined by inequality constraints was equivalent 
to the saddle value of the Lagrangian function of that pro- 
gramming problem. The interpretation of the Kuhn-Tucker 
results when applied to economic theory was consistent 
with concurrent economic research efforts. 

In 1961, R. Pfouts [7] applied the Kuhn-Tucker neces- 
sary conditions to a convex problem formulation of cost and 
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production in the multi-product firm in a perfectly com- 
petitive economy. In this formulation, the multi -product 
firm was modeled as a homogeneous economic entity rather 
than an aggregation of single product firms, as in the 
Hicksian model. This approach characterized each pro- 
duct as competing for a portion of the total amount of 
fixed factors of production (long run sense) assessing 
a "set-up" or transfer cost for changes in product mix 
which required corresponding changes in the allocation 
of the fixed factors between periods (short run) . Levels 
of the variable factors of production (short run sense) 
were assumed to change with the output level. 

The model seeks to minimize a total cost function 
consisting of variable, fixed, and transfer cost terms 
constrained by technology and continuity of fixed factor 
allocation to specific products. The convex problem is 
solved by applying the Kuhn-Tucker necessary conditions 
yielding the following results at equilibrium. Consis- 
tent with the Hicksian model, the imputed value of mar- 
ginal output must be at least as great as the marginal 
costs of that output. But additionally, the imputed value 
of the reallocation of fixed resources from one product 
to another must at least be equal to the value of the in- 
creased marginal product resulting from the change minus 
the cost associated with the marginal transfer of those 
fixed resources. Thus, only when an excess of all fixed 
factors exists can the multi -product firm be considered 
an aggregation of single product firms. 
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Four years later, T. Naylor [7] expanded Pfouts' anal- 
ysis of the multi-product firm in a competitive economy to 
a monopolistic-monopsonistic economy. Maintaining Pfouts' 
concept of a multi-product model, Naylor optimizes the 
economic dual of that model. At equilibrium the decision 
rules characterizing the optimum substantiated Pfouts' 
previous conclusions. 

In 1958, G. Dantzig and P. Wolfe [9] introduced a 
decomposition principle for linear programs. The princi- 
ple is applicable to programs which have some constraints 
which "bind" together various sets of constraints which 
are independent of each other. The original problem is 
decomposable into subprograms, a characterization of each 
muepeiideut tens uiaiut , and a master program, an aggre- 
gation of the subprograms. A full master program is formed 
from the convex combinations of all the extreme points of 
the feasible regions of each of the independent constraints. 
Each basic feasible solution to the master problem generates 
dual variables. A function of the dual variables establishes 
an optimality criterion for the value of the master program 
objective function. If this criterion is not met for any 
one or more of the independent sub^roblems, the subprogram 
for that constraint set is optimized providing a new vec- 
tor to enter the basis and generating a new basic feasible 
solution to the master program. This procedure continues 
iteratively until the optimality criterion is achieved for 
all independent subprograms or it is evident an optimal 
solution does not exist. 
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In 1958, L. Hurwicz [10] describes the "greed process" 
as a model of information decentralization for a class of 
decomposable economic environments. In this context de- 
composable means free of external (dis) economies of scale. 
The "greed process" embodies the intuitive concept of de- 
composition because each economic agent is assumed to be 
concerned with his actions alone, even though this may be 
against his economic best interests in the operation of the 
system. Each agent lacks any knowledge of the internal 
functioning of other agents in the economy. Hurwicz proves 
that the "greed process", the interaction of these inde- 
pendent economic agents, achieves Pareto optimal conditions 
but requires more information to achieve these conditions 
than the neoclassical competitive market mechanism. Al- 
though this "greed process" is shown to be applicable to a 
broader class of economic environments, it lacks the infor- 
mational efficiency of the perfectly competitive mechanism 
which decentralizes the economy by prices alone. 

In 1963, Y. Ijiri [11] applied a special form of lin- 
ear programming, "goal programming", to accounting models 
of the firm. Goals may take the form of lower bounds on 
production or upper bounds on input usage. The procedure 
emphasizes the decomposition of definitive management goals, 
which correspond directly to profit, into a set of subgoals 
which are more tractable for subordinate agencies within 
the firm. Management seeks to minimize the firms deviation 
from the firm's goals. The deviations from the firm's goals 
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are an aggregate of the subordinate agencies' deviation 
from their respective independent goals. If subordinate 
goals are compatible, a case in which all goals can be 
satisfied simultaneously, specific goal levels and devi- 
ations from those levels are sufficient information to 
decompose the firm and iteratively achieve an optimum 
solution. If subordinate goals are incompatible, a case 
in which all goals cannot be satisfied simultaneously 
but some goals must be satisfied sequentially, Ijiri 
proposes a "preemptive priority", an ordering of goal 
achievement. Further, if there are several goals of the 
same order, a weighting scheme can be used to provide 
relative differentiation among goals of equal priority. 
Utilizing a procedure similar to the simplex algorithm 
of linear programming, the goals with the "lowest cost", 
i.e., highest priority and weight, sequentially enter the 
basis until an optimum solution for the existing priority 
and weighting scheme is achieved. 

Recently, T. Ruefli [12] analyzed resource allocation 
in a three level firm by utilizing a generalized goal de- 
composition model. The organization consists of a cen- 
tral agency, management and operational units. Each level 
is vertically interrelated by a specific information flow. 
Levels of organization are assumed to be horizontally in- 
dependent. The central agency seeks to maximize the im- 
puted value of scarce resources subject to the allocation 
of those resources to different managements. Managements 
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seek to minimize the deviations from these specified goals 
subject to the productional relationship of their indepen- 
dent operational units. Managements pass the shadow prices 
of their goals to subordinate operating units and the cen- 
tral agency. Each operating unit minimizes it's management's 
imputed value of their (each operating unit's) production 
activity subject to each operating unit's particular tech- 
nology. Operating units propose activity levels at that 
shadow price to their respective managements. Concurrently, 
the central agency determines revised goals. Managements 
sequentially determine revised shadow prices given these 
new goals and new production activities. The procedure 
continues iteratively until a desired optimum is achieved. 

In 1953, C. Carrol [13] proposed the Created Response 
Surface Technique, a penalty function approach to the solu- 
tion of nonlinear constrained mathematical programs which 
modeled complex industrial processes. The technique 
"created" a form of weighted penalty function from the 
original objective function and original constraints. The 
procedure seeks to optimize the unconstrained created func- 
tion in a sequence of steepest ascent steps. Each itera- 
tion of the procedure seeks the optimum to an artifical 
objective surface. This optimum in turn generates another 
artificial objective surface whose optimum value is sought. 

As the number of iterations increase beyond all bound, the 
artifical objective surface approaches the surface of the 
original objective function and the artifical optimum ap- 
proaches the original optimum. Throughout the procedure 
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feasibility is assured by a weighted "boundary-repulsion" 
term which is sequentially reduced at each iteration. 

In 1961, P. Wolfe [14] formulated a dual program for 
a nonlinear, convex, differentiable, primal, objective 
function minimized over a convex constraint region. 
Applying the Kuhn-Tucker Equivalence Theorem conditions 
to the Lagrangian Function of the primal problem, Wolfe 
proved the existence of the dual solution at optimality. 
Also at optimality the values of the primal and dual ob- 
jective functions are equal. 

In 1968, M. Balnski and W. Baumol [15] interpreted 
Wolfe’s dual formulation to a nonlinear primal problem of 
profit maximization constrained by technology in a com- 
petitive economy. Maximization of a concave profit func- 
tion subject to concave technological constraints resulted 
in the dual minimization of a Lagrangian function consist- 
ing of a fixed profit function, the value of unused re- 
sources and the total marginal opportunity losses of all 
outputs. The minimization of the dual objective function 
minimized the value of unused resources and the total 
marginal opportunity losses of all outputs. Regrouping 
of terms in the dual objective function provides a mathe- 
matical interpretation of economic rent which is the dif- 
ference between total profit and the imputed value of all 
inputs. The existence of economic rent is attributed to 

the concave object function representing diminishing re- 
* 

turns. This assumption indicates that the imputed value 
of inputs should be less than total profit. 
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In 1963, A. Fiacco and G. McCormick [16, 17, 18] 
developed a modification of Carroll's Created Response 
Surface Technique which solves a constrained, nonlinear 

i 

mathematical program in a sequence of unconstrained mini- 
mization problems. The original programming problem is 
to minimize a convex objective function subject to con- 
cave constraints. A penalty function is formed from the 
original objective function and constraints. The sequen- 
tial unconstrained minimization of this form of penalty 
function approaches the optimal value of the original 
problem by iterative gradient search methods. A "boun- 
dary repulsion" term insures feasibility from an interior 
starting value. Primal-dual feasibility, recommendations 
b cis 6 d czp c n n o^p db d ti i on 3.1 sxp science cidd tibc thcorsticcl 
structure of the technique are presented. 
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Ill . THE CREATED RESPONSE SURFACE TECHNIQUE 



A. PROCEDURE DESCRIPTION 

In activity analysis, mathematical programs are uti- 
lized to model the economic activity of the firm. Cus- 
tomarily, the model consists of an objective function, a 
function of endogenous variables under the control of 
management, and constraints, levels of exogenous variables 
which are determined from the economic environment. As- 
suming neoclassical rational behavior, the objective func- 
tion is maximized or minimized, as appropriate within the 
bounds defined by the constraints. 

If the objective function and constraint equations 
are linear functions, linear programming techniques may 
be used to determine the appropriate decision rules or 
levels of endogenous variables. If the objective func- 
tion and constraint equations are nonlinear but separable 
functions, separable programming techniques may be used to 
solve the program. Similarly, quadratic programming 
methods may be used if the objective function is quadratic 
and the constraint equations are linear. 

There are few developed techniques of solution for 
those models which are not classified as linear, quadratic, 
or separable. However, an auxiliary function method has 
been developed for this class of models. The general ob- 
jective of this approach is to transform the original con- 
strained problem into an unconstrained "auxiliary function" 
which may be optimized by several available techniques. The 
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auxiliary function incorporates the objective function and 
constraint equations of the original problem. 

One particular class of auxiliary functions consists 
of the original objective function and a penalty term 
which algebraically contributes a penalty for violation 
of the constraints of the original problem. This penalty 
term may be logarithmitic or reciprocal -linear depending 
on the particular model. The intent of this penalty func- 
tion approach is to insure that the number of infeasible 
solutions is decreased as the procedure approaches the 
optimum of the auxiliary function. The optimum of the 
auxiliary function will be the optimum of the original 
objective function if the procedure is convergent. 

The Created Response Surface Technique [13] is one 
form of the penalty function approach to the auxiliary 
solution of a constrained mathematical program. 

As an explanation of this technique, consider the 
following problem: 

MAX: f (x-^ • • *Xj • • *Xj) 

(A) 

S.T. gi (x) > k. i = 1,. . . ,1 

where , 

f(x) is a concave function with continuous first and 
second partial derivatives, 
g^ (x) for i = 1,...,I are concave functions with con- 
tinuous first and second partial derivatives. 
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The constraint region formed by the (x) for i = 1,...,I 
is a convex region. If the nonnegativity of the dependent 
variables is required, these would be included as addition- 



constraint equations, ( i = I+l 


)•••>! > 


in the problem 


The 


created response function 


problem for (A) is: 


(B) 


MAX: P(x,r) = f(x) - r 


I 

y u i 






i=i Si 00 


- k. 
i 


elements of this function are: 






1. 


Created response function 


: P(x,r) 




2. 


Objective function: 


f(x) 




3. 


i^k constraint equation: 


gi(x) > 


k. 

i 


4. 


i^h constraint level: 


k. 

1 




5. 


, t ll -v. n n 4 w ^ n 1 4 A..4 a + 4 A A . 

-L i UCXpAUtax Uw V XU4.0-VU . 


1 fry C v ^ 


_ V 


6. 


i^* 1 subjective weighting 
factor : 


w - ; w. 

i * i 


> 0 


7. 


Penalty term: 


i 

y 


w . 
i 
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(x) - k. 


8. 


Penalty term weighting 
factor: 


r ; r > 


0. 



The original constrained maximization problem has been 
transformed to an unconstrained maximization problem. Sev- 
eral methods are available to perform this unconstrained 
maximization. The second order gradient search method will 
be used because of its proven computational efficiency. 4 

4 

Fiacco, A. and McCormick, G., "Computation Algorithm 
For the Sequential Unconstrained Minimization Technique for 
Nonlinear Programming," Management Science, Vol. 10, No. 4, 
p. 607, July 1964. 
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Intuitively, the maximization of the created response 
function should decrease the value of the weighted penalty 
term relative to the value of the objective function. The 
procedure seeks to maximize the function and avoid the 
bounds of the constraint region simultaneously. If any one 
constraint approached its binding value, g^(x) = k^, the 
penalty term would increase beyond bound. As the weighting 
factor of the penalty term, r, is sequentially reduced to 
zero, the procedure will allow constraints to become "more 
binding" in the sense of approaching their respective boun- 
dary values. In the limit, the value of the optimum of the 
created response function approaches the value of the opti- 
mum to the original problem. 

Assume that the I subjective weighting factors have 
been assigned. Starting the procedure with an initial fea- 
sible point, x , and initial value of r, r p a response 
surface, P(x,rj), is created. The second order gradient 
search method seeks the optimum value of that response 
surface starting at x r The optimum value of PCx,^) 
yields a new starting point, x 2 . P(x,r 2 ), where r 2 < r x , 

is a new created response surface whose optimum value, x 3 , 
is found by using second order search starting at x 2 . Then 
r 3 is chosen with r 3 < r 2 and the process continues. Each 
new r value creates a new response surface until as, 

Lim r c ■+ 0 
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the optimal value of the £ th created response surface ap- 
proaches the optimum value of the original problem. 5 

B. PRIMAL FORMULATION 

1 . Practical Application 

The Created Response Surface Technique exhibits 
several characteristic properties which indicate the fea- 
sibility of practical applications of the technique in 
activity analysis. The following description of these 
properties will be used to sketch the procedure's poten- 
tial application and to state these properties without the 
rigorous proofs developed in References [16, 17, 18]. 

The Created Response Surface Technique could be 
used to solve models of complex production processes with- 
out requiring that linear, quadratic, or separable con- 
ditions be met. Linear models assumed constant returns to 
scale. Nonlinear models could permit production functions 
which demonstrate increasing, constant, or decreasing re- 
turns to scale over the range of the endogenous variables. 
Linear models assumed production processes were independent 
and in fixed proportion; nonlinear models could introduce 
the interactive effects of different production processes 
at various levels of mix which could not be described by 
quadratic or separable functions. 



5 

Fiacco, A. and McCormick, G., "The Sequential Uncon- 
strained Minimization Technique for Nonlinear Programming, 

A Primal-Dual Method," Management Science , Vol. 10, No. 2, 
p. 361, January 1964. 
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The convergence of the optimal values of (A) and 
(B) is assured in the limit. The stability, which is the 
local improvement of the created response function at each 
iteration, provides assured improvement of the objective 
function for a finite number of iterations of the procedure. 
A criterion for local improvement from an initial starting 
point may be achieved. 

Constraint feasibility is assured in all itera- 
tions of the procedure. Although this characteristic may 
require an increased number of iterations to reach a solu- 
tion which may be achieved in less iterations by permitting 
temporary infeasibility, the temporary infeasibility may 
be economically uninterpretable in the context of the orig- 
inal problem. 

By utilizing gradient search techniques to opti- 
mize the created response function, local improvement for 
each value of r is achieved. Assuming that the response 
surface is well behaved, changes in the values of the en- 
dogenous variables should reflect local changes during each 
iteration. These endogenous variables may represent phy- 
sical production processes. Minor changes in these pro- 
cesses represent "local" improvemeht which may be less 
costly to introduce them than massive changes required if 
the global optimum is achieved indirectly without inter- 
mediate steps. The intermediate "optimum" values of the 
variables generated by each iteration of the search tech- 
nique may facilitate making physical changes in the appro- 
priate production process. 
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